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ABSTRACT 

We consider the eiectromagnetic response of a topoiogicai Weyi semimetai (TWS) with a pair of Weyi nodes in the buik and 
corresponding Fermi arcs in the surface Briiiouin zone. We compute the frequency-dependent compiex conductivities (cb) 
and aiso take into account the modification of Maxweli equations by the topoiogicai S-term to obtain the Kerr and Faraday 
rotations in a variety of geometries. For TWS fiims thinner than the waveiength, the Kerr and Faraday rotations, determined by 
the separation between Weyi nodes, are significantiy iargerthan in topoiogicai insuiators. In thicker films, the Kerr and Faraday 
angles can be enhanced by choice of film thickness and substrate refractive index. We show that, for radiation incident on a 
surface with Fermi arcs, there is no Kerr or Faraday rotation but the electric field develops a longitudinal component inside the 
TWS, and there is linear dichroism signal. Our results have implications for probing the TWS phase in various experimental 
systems. 


Introduction 

In recent years, condensed matter physics has witnessed the emergence of novel quantum phases characterized by topology 
rather than by symmetry breaking. The best studied of these are the topological insulators, which have a bulk band gap, but 
with gapless edge or surface states protected by time reversal symmetry and characterized by topological invariants. More 
recent predictions suggest that nontrivial topological properties can also arise in certain systems whose gapless band structures 
are characterized by point or line nodes.A particularly interesting state of matter is the topological Weyi semimetal 
(TWS).^’®’^ These are phases with broken time-reversal or inversion symmetry, whose electronic structure consists of pairs of 
Weyi nodes, points in the bulk Briiiouin zone (BZ), which are at the chemical potential and act as sources and sinks of Berry 
curvature. This is predicted^ to lead to unusual surface states that are gapless on disconnected Fermi arcs with end points at 
the projections of the bulk nodes onto the surface BZ. 

A number of candidate material systems that should exhibit a TWS phase have been proposed. These include pyrochlore 
iridates A 2 lr 207 ,^ spinels® and multilayers of topological insulators and trivial insulators,^** all of which break time rever¬ 
sal. There is a recent report of the experimental observation" of surface Fermi arcs in TaAs, which is a TWS by breaking 
spatial inversion. In addition there are many theoretical predictions about the unusual transport and magnetic properties of 
TWSs 8,10,12,13,13-19 

In this paper, we theoretically address the electrodynamic response of a TWS with broken time reversal symmetry, fo¬ 
cusing on Kerr and Faraday rotations and linear dichroism. We show that our predictions are sensitive to four nontrivial 
characteristics of a TWS: (i) the topological Weyi nodes which lead to nontrivial a^y in the absence of an applied field, (ii) 
the nodal excitations leading to optical conductivity C7xr(®) ~ 0), (iii) the unusual surface states with Fermi arcs, and (iv) the 
modification of Maxwell equations inside the TWS via a topological 0-term. 

In brief, our results are as follows. We find that there is a very important difference between the EM responses for radiation 
normally incident on (A) a surface that does not support Fermi-arc electronic states, versus (B) a surface that does. Specifically, 
for a pair of nodes separated along the k^-direction in the bulk BZ, the (x,y)-plane has no Fermi arc states; see Fig. 1. We 
find Kerr and Faraday rotations only in case (A), i.e., for light incident on the (.r,y)-plane. For case (B), there is no Kerr and 
Faraday rotations, but the electric field develops a longitudinal component inside the TWS and there is a finite magnetic linear 
dichroism signal. 

We look at various experimental geometries - a TWS film thinner than the wavelength of light, a semi-infinite slab and a 
film with thickness comparable to wavelength - and identify cases where large measurable signals are obtained. We discuss at 
the end of the paper, various experimental systems where our results can be tested. In addition to the TWS materials already 
mentioned above, our results may also be generalized to recently discovered Dirac semimetals, where a magnetic field B 
separates the Weyi nodes. 





Figure 1 . (Color online) (a) Schematic k-space picture showing the location of Weyl nodes (blue and red dots with outward 
arrows as sources and inward arrows as sinks of Berry curvature) along axis in bulk Brillouin zone (BZ) and the Fermi 
arcs (lines ending at the projection of the Weyl nodes) on the surface BZ’s, shown as blue squares. Note that surfaces 
perpendicular to z axis have no arcs, (b, c) Real-space geometry of slab and electromagnetic waves. The blue (dark) surfaces 
of slabs support the arc surface-states and yellow (light) surfaces are ones without surface states. In (b) we show the Kerr 9k 
and Faraday Op angles. For geometry shown in (c) there are no Kerr and Faraday rotations, instead E acquires a longitudinal 
component inside the Weyl semimetal. 


Results 

The conductivity tensor 

The first step in computing the Faraday and Kerr effects for a TWS is to obtain the conductivity tensor Gap (w) = (7^^ + 
using the Kubo formalism. To make the physical ideas clear, we focus on a TWS with two Weyl nodes with opposite chiralities, 
located at ±b = ±(0,0,/?). The low-energy physics is described by the linearized Hamiltonian 

//(k) = ±riyfa-(k±b), (1) 

where k < kc, the momentum cut-off, vp is the Fermi velocity, and a = (Ot, O)., ov) are Pauli matrices. 

In addition to the known result^^ = e^\co\/6hvp, we also determine While the linearized H is sufficient 

to obtain g'{co) for o <C Wc = vpkc, it leads to a pathology in a”{co), even in the low frequency regime. This is avoided by 
using a lattice regularization. We find that has a spurious — 1 /o) term for the linearized dispersion, which is in fact exactly 
cancelled by the diamagnetic term in the Kubo formula. Such a diamagnetic term is indeed present for the lattice H, but absent 
for a strictly linear dispersion; see Methods for more details. 

Using e = e/Eo = £b + i<7xt/eoW, where Eb is the bound charge contribution, we then find the real and imaginary parts 
of the dielectric function. These are given by e'{co) = Eb + (ac/37rvf )ln |(a)^ — 4a)^)/a)^| and e" = acj^nvp, where a = 
/AtzeoHc is the fine structure constant. 

For the Hall response, we find the aP' correction to the known d.c. result^’ and obtain <7^(w) = e^b/Kh+ [e^b/6nh)(>P/ — 

v^b^)- The gP terms in e' and are shown here for the linear model; for the lattice model the prefactors, which must be 
found numerically, have essentially the same structure with (Oc ^ vp/a. We note that disorder produces subleading correc¬ 
tions^^ to (7^, and does not^® affect . Therefore the above results remain valid at least for weak disorder. 

Electromagnetic Response 

We now consider in turn the electromagnetic response of a TWS with light incident on: (A) a surface that does not support 
Fermi-arc states, which is the (;ic,y)-plane with nodes separated along the z-direction, and (B) surfaces that do support Fermi- 
arc states; see Fig. 1. We look at three geometries: (1) In an ultra-thin film, which is thinner than the wavelength, the a^p (co) 
of TWS impacts the boundary condition at the interface between two topologically trivial media. (2) When the TWS is a 
semi-infinite slab, we need to consider the modification in the Maxwell equations arising from the topological 0-term. (3) In 
a film of finite thickness comparable to the wavelength of light, we take into account both the modified Maxwell equations 
and interference phenomena arising from multiple reflections. 

Kerr and Faraday rotations 

• Case (Al): We first consider light incident on the (.r,y) surface with no arc states [see Fig. 1(b)] for an ultra-thin film 
with thickness d with a <C t/ <C A, the wavelength of light. In this limit, the nontrivial properties of the TWS enter only 
through the boundary condition for e.m. fields in the two non-topological media on either side (z < 0 and z > 0). The 
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surface current density at z = 0 is with the surface conductivity tensor O'^p — dOap for the TWS thin film. 

The rotation of the polarization in both reflection and transmission is governed by <7^. 

Using right and left circularly polarized transmitted fields =E^^zL iE\, the Faraday rotation is given by Qp = (arg£''|_ — 
argiiL )/2. An analogous treatment of the reflected fields yields the Kerr rotation. We thus And 


tan0f = 2abd/n\2 +a(od/3vp\ 

(2) 

^ —Aabd/n 

^ {a(od/3vpy+ 2acod/3vp-f{2abd/ny' 

(3) 


Here we ignore the o-dependence of and show results for free standing Aims. 

At the lowest frequencies, we And tan Qp K,bda/7i = — cot Qp- Note that Op for a TWS is enhanced by a factor bd^ 1, 
relative to the small (< 1°) Faraday rotation for the surface of a topological insulator with time reversal breaking. 

We And that Op for a TWS can also be very large. For instance for typical Aim thickness of d = 5Q —100 nm it reaches 

values of 70—80°. 

We also note that there is a flnite frequency regime where the results are independent of Aim thickness d. For 6bvp/% < 
CO < 3vp/ad, two second-order terms in denominator in eq. (3) can be ignored. We thus obtain tanO/r = —Evpbjnco, 
independent of d. For a typical values of d = 100 nm, vp = 10® m/s, the corresponding photon energies should be 
1 < hco < 5 eV, corresponding to the visible range, which is accessible in experiments. 

• Case (A2): For a semi-inflnite TWS (z > 0) with light incident on a surface without arc states [see Fig. 1(b)], we must 
take into account the modiflcation of the Maxwell equations inside the TWS. The topological axion term^^’^^ in the 
action with 0 (r) = b • r leads to 

V-D = p + K-b-B, (4) 

VxH = dD/dt+J-KhxE (5) 

with K= {2al7z)y^eoliM)=e^l7zh. 

In our geometry [Fig. 1(b)] the E and B flelds are in the (x,y) plane. As usual, B = pH with jj. = p./jjp, = 1 + Xm- The 
effective dielectric tensor (a, j3 = x,y) of the TWS can be written as 

£a,p = £Sa,p - {2aln)icblco) ( t >)„ ^ . (6) 

£ is the complex dielecti'ic function £ = eh + iOxx/soCO and the Pauli matrix term arises from b x E w'lihh — bz. Thus 
the constitutive relation for a TWS is of the gyrotropic form^®’^® D = £E — /g x E, well known in the electrodynamics 
of ferromagnetic materials, with the nodal separation b as the gyration vector”® g. 

We use the Jones^® basis of eigenvectors of to obtain the complex refractive indices 


= \ e±aXob/ii^ 


(7) 


for left and right circular polarized light inside the TWS. with Aq = 2ncl(O the wavelength in vacuum. The difference 
between leads to birefrigence. 

The tangential electric flelds are continuous across the vacuum-TWS interface, because the total surface current density 
at z = 0 vanishes. First, there is no transverse response localized on the z = 0 surface. Second, in contrast to metals, the 
longitudinal current density in a TWS is not localized at the surface, given the long penetration depth 5 ~ e'X / ne" ^ A 
over which flelds decay in the TWS. 

Writing (n^ — «+)/(!— n^rT ) = r\e"^, we And the Anal result for the Kerr rotation to be (see Methods) 


Qk 


1 _[ /2T7sin^\ a^c Xb 

2 '’^” ) ~ ^Tt^vp e'{e' -ly 


( 8 ) 


where A = Ao/v/U is the wavelength inside the TWS. 

The dielectric constant of TWS can be rather large, e.g., in a topological insulator (Tl)-ordinary insulator multilayer 
Bi 2 Se 3 has £' ~ 30— 80. Therefore the Kerr rotation from a single interface is small ~ 10^® — 10^^ rad (like in 
Ti’s 21,22,27,28) show next that reflection from a thick Aim can substantially enhance the Kerr rotation. 
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• Case (A3): We next consider a TWS slab of thickness d and dielectric function £2 + sandwiched between two non- 
topological media with refractive indices ni = 1 (vacuum) and na (substrate). As before, we look at circularly polarized 
light incident on a TWS surface with no arcs [Fig. 1(b)]. Using eq. (7) and na = refractive index for the TWS 

is given by — n 2 (l where = ±aAo^/27r^£2 + *£2 

The multiple scattering from the two interfaces lead to interference effects with reflection and transmission coefficients 
that vary with thickness d. Here we mention only some of the results focusing on the cases where a large Kerr/Faraday 
rotation is predicted. At reflection maxima, where d = pX jl for an integer p, the Kerr rotation for a free standing slab 
with n\ = ^3 is given by 


Ok — — tan ^ 




(9) 


independent of d. For THz radiation Ao = 10/im and Fermi velocity vp — we get large dp — —it 12. We can 

also get a large Kerr angle ^njl — n^d^ / aX^jb near a reflection minima, when d = (2p + \)X/A, provided one can 
also satisfy n 2 — by choice of substrate. The Faraday rotation of transmitted waves can be similarly enhanced 

for refractive indices satisfying 112 — ^Jnpvi, with 


Bf 


n ^ aXfjb 
2 d,nn2 ’ 


( 10 ) 


where the second term is very small for wavelength of about Xq < 10/rm. 


Linear dichroism 

Finally we turn to light incident on a surface of the TWS that supports Fermi arcs. We take the incident wave propagating 
along X [see Fig. 1(c)], the Weyl nodes are at ±hz in k-space with Fermi arcs in the {ky,k^) surface Brillouin zone. We will 
focus here on Kerr reflection from a semi-infinite TWS (.r > 0). 

Consider the propagation of a plane wave with wavevector k = nco/c inside the TWS. The modified Maxwell equations 
(5) lead to 



where £ = £* + i<7xx/£o(0 and e'^ = 2acbjn(0t. In order to have a solution, the refractive index n must be either (i) = £ 

or (ii) n\ = e — d^je. As we shall see next, the subscript here refers to whether the incident E held is || or _L to the node 
separation along z. For case (i), E = and matching boundary conditions at x = 0, we And no rotation of E 

upon entering the TWS. This is easy to see since b • B and b x E in eq. (4) are identically zero, when E| |z and H| |y. 

Case (ii), relevant for incident E||y and H||z [see Fig. 1(c)], is more interesting. Now the solution inside the TWS is 
E=(y — /x d„/d) 9. Thus we And no Faraday or Kerr rotation of the polarization in the (y,z)-plane, however, 

we do see that E acquires a longitudinal component inside the TWS. These conclusions, which are further elucidated below, 
can be reached just by looking at the solution of eq. (4) inside the TWS (x>0). The boundary conditions at x = 0 (discussed 
in Methods) only determine the magnitude of Eq. 

First, we comment on the absence of Kerr and Faraday rotations in our geometry. This is related to the fact that the off- 
diagonal £yj in eq. (11), which could have caused rotations in the (y,z) plane, vanishes. Although time reversal symmetry is 
broken, there is no yz Hall response. In a very different set up, where an external E B term creates a charge imbalance between 
nodes, one can get non-zero Ey^. This leads to a Faraday rotation,which is independent of node separation and vanishes in 
the limit of equal Fermi energies at the two nodes. 

Next, we remark that the longitudinal component of E obtained above is, in fact, just the consequence of broken time 
reversal (and exists even in magnets^®); it is unrelated to topological properties of the TWS. The latter can be probed in 
magnetic linear dichroism (MLD) measurements on surfaces with arcs. In cases (i) and (ii) discussed above, the waves 
propagate in the TWS with different refractive indices n|| and np. MLD measures the difference between absorption of light 
polarized linearly in different directions.^® Thus the MLD signal yields 


n 


n 

X 



/2a \ ^ / bX'\^ vp 

\ n ) \2n^) c 


( 12 ) 
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Discussion 


We have presented in this paper detailed predictions for the Kerr and Faraday rotations and linear dichroism in a Weyl 
semimetal. Our results give insights into various aspects of the TWS state including its topological characteristics; Kerr 
and Faraday rotations from surfaces with no arcs and linear dichroism signal from surfaces supporting arcs. Perhaps the 
simplest candidate for testing our predictions is a topological insulator-trivial insulator multilayer geometry,'** which has only 
two Weyl nodes. For the TWS phase in the complex materials one must take into account multiple pairs of Weyl nodes. This 
can be done by a linear superposition of the results for individual pairs. Candidate materials for exploring the TWS phase 
include pyrochlore iridates A 2 lr 207 in bulk^’^® and in [111] thin films,and HgCr 2 Se 4 spinels. 

We should also note other systems whose ground state is not a TWS, but where Weyl nodes can appear once time reversal 
is broken by an external magnetic field. In the recently discovered Dirac semimetals Cd 3 As 2 ^^“^^ and NaaBi,^®"^^ each Dirac 
node is expected to split into two Weyl nodes with a separation that grows linearly with B. On the other hand, in systems with 
a quadratic band touching,^® the application of a magnetic field leads to Weyl nodes with a separation proportional to \/B. Our 
results on Kerr, Faraday and linear dichroism signals can also be generalized to these systems. 


Methods 

We present here detailed derivations of many of the results reported in main text. We begin with the Kubo approach to the 
longitudinal and transverse conductivities of a TWS. We then turn to transmission and reflection matrices connecting scattered 
waves to incident waves. In the next two subsections, we discuss the modified Maxwell equations for a Weyl semimetal (TWS) 
and reflection from a semi-infinite system. Finally, we give the details of the calculations for the case when light propagates 
perpendicular to node separation. 

We use following notation in the text and Appendices. For all quantities X — X' + iX" with X' — ReX and X" = ImX. 
Subscripts 1,2 etc. refer to medium 1,2 and not to the real/imaginary parts. The permittivity e and the permeability p. are 
defined via the constitutive relations D = eE and B = /IH. Their “relative” counterparts, which are dimensionless, are defined 
by e = e/eo and p — p/po, where eo{po) is permittivity (permeability) of vacuum, e is the dielectric constant which then 
defines the refractive index as n = ^/e, and p is connected to magnetic susceptibility Xm p = I + Xm- We assume that the 
susceptibility of the TWS is negligible. 


Optical conductivity tensor of TWS 

In this subsection we calculate the longitudinal Oxx and transverse Oxy conductivities of a TWS. We consider the simplest case 
with only two nodes located at ±b, where we choose b = bz. Near the nodes, the linearized Hamiltonian is given by 

H{k) = ±hvFd-{li±h). (13) 


The symbols used are defined in the main text. We obtain the conductivity from the Kubo formula 

(yal3i(0) = -^limTl^l3iq,C0). (14) 

CO q->0 ^ 

Note that there is no diamagnetic term for a strictly linear dispersion, a point which we will come back to below. The 
polarization function (q, w) is given by the current-current correlation function 

na/3(qF'®«) = dTe‘“'’\Trfa(.q,T:)fpi-q,0)) (15) 

where 'f' is the volume of the system and the current density operator 
5H 

J = -^=±evF(J. (16) 


The real frequency behavior is obtained by analytic continuation i(0„ —W -I- /0+. For each node we obtain 

/(e,.(k))-/(e,-(k)) 


na,/t(®) = 


2 2 
e 


y 


^ y _ 

hco + £;/ (k) - e;(k) -f /0+ 


(k/|(7a|k/')(k/'|(7j3|ki), 


where f{x) = l/(e^^-)-1) is the Fermi function. The quasiparticle energies and eigenstates are obtained from //(k)|k/) = 
e,(k)|k/) with i= 1,2 labeling two eigenstates at each wavevector. 
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We evaluate the longitudinal and transverse polarizations (o) = 11^^ (o) + (o) when the Fermi energy lies with 

nodes. Therefore, at the low frequency limit (O < (Oc the total longitudinal conductivity from both nodes is given by; 


, CO 

OTT/Z Vf 

= - 


3n^hvp 


CO 


CO , 


-In 


co^ - 4®/ 


CO'^ 


(17) 


(18) 


Note that the spurious 1 /® term in <7^ is an artifact of effective linear Hamiltonian with no diamagnetic response. We 
have shown that this problem is cured when we consider a lattice Hamiltonian®’^® 


//(k) = 2f sin(^va)(7i + 2f sin(A:ya)(7v + 2f [2 + cos(hfl)](7j — 2f [cos(^ja) +cos{kya) +cos(^ja)] ffj, 


where a is the lattice spacing. For the lattice model we found that the diverging 1 /co terms in paramagnetic and diamagnetic 
responses precisely cancel each other. Therefore, in the properly regularized model, we are left only with log term in the 
equation above for cy^{co) with a cutoff that is effectively cOc t ^ hvp/a. This lead to a log term in the real part of dielectric 
constant e'(®). 

In the limit of low frequencies, we hnd the transverse conductivity is given by: 




CO 


(19) 


where we denote here by kc, the cut-off along k^ axis. 


Transmission and refiection coefficients for uitra-thin fiims 

In this subsection we derive the reflection and transmission matrices in the ultra-thin him limit where the TWS him thickness 
<i <C A, the wavelength of light, even though d ^ a, the lattice constant. Thus the two media i {z < 0) and j (z > 0) are 
separated by a TWS which gives rise to the boundary conditions at z = 0. 

Taking the incident wave to propagate along the z axis in medium i as E"' = {Eljc+Eyy)e‘^‘^', the expressions for EM waves 
in the two media are as 


Medium i: E,- = 


4 


+ e 


-ikjZ 


Medium j: E,- = 




( 20 ) 


where wavevector is k,- = /c using notation n, = y/ei, and superscripts r and t stand for rehection and transmission compo¬ 
nents. We assume medium i is vacuum n, = 1 and medium j (substrate) has refractive index rij = n. In the main text we only 
quote the results for a free-standing him, i.e., for tij = \. 

The EM waves should match on the boundary and are given by n x (E; — Ey) = 0 and n x (H, — Hj) = Jj, where n = —z 
is normal to interface. The surface current density Js = (7®E, where the surface conductivity tensor = da^p with the bulk 
a^p for the TWS given by the results in Appendix . 

Solving the boundary value problem, we get the following expressions for transmitted wave 


Ei = 


2(1 +n + -^^a^) 


(1 -fn H—1-(J* )^ -f ( —(J® )^ 

V® ^ c£o -tit ^ 


C£q 


Ei = 


— O’® 

ceo 


fl-l-n-l-1-C7® )^-|-( —(7® )^ 

1 -r « -r g ^ \ cEn 


7 El. 


( 21 ) 

( 22 ) 


We use these results to calculate Earaday rotation in Eq. (2) in main text. Similarly, for rehection components we have 
1 — fn-I--^(7®'j — f-l-C7l') 

r ^ V xxj ycEp xyj 

^ (l+„ + J_0®)2 + (J-o®)2 

— O® 

_ _ eEq ^ 

“ (l+« + ^i^C74)2 + (^i^C74)2 


(23) 

(24) 


These relations are used to calculate the Kerr rotation in Eq. (3). 
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Modification of Maxweii equations in TWS 

The low energy electromagnetic response of Weyl semimetal is described by an spatially varying axion term. Including 
conventional Maxwell term, the full action S = So + Sg of system is^' 

So = J d\dt I - -^PpivP^'' + \ppiv I > 

5e = ^/r/^x*{0(r,f)e^''“^iVv^’a/t}. (25) 

Here the tensor stand for electric polarization and magnetization as = cP‘ and The axion field 0 

varies with space and time as 9{r,t) = 2b • r — 2bot, where b(fio) denotes separation of nodes in momentum (energy) space.^^ 
We set bo to be zero, since in the problem of our interest both the Weyl nodes are at the same chemical potential. 

Varying action with respect to A^, the equation of motions are obtained as follows. 

- + -^e^^^PdyiOF^s) = . (26) 

Mo 27ZHo 

These equation of motions yield the modified Maxwell equations in the presence of axion field in the effective action.In a 
more conventional form, they are as 


VxE + ^ = 0, 

dt 

VB = 0, 


V D = p 


2a 


-b B, 


VxH = 




Mo 

2a 

71 


-bxE, 


(27) 

(28) 

(29) 

(30) 


In the absence of topological terms proportional to b, we assume that the constitutive relations are isotropic with D = eE 
and B = /iH. Here e — e/eo = £b + iOi/eoCO and ji — jl/^ = 1+ Xm, where £), denotes the bound charge contribution and Oi 
is longitudinal conductivity. With the topological term, the dielectric constant acquires off-diagonal elements, and is given by 


V = el-— 

K CO 


where t'" is Pauli matrix that comes from writing b x E in terms of its components when b = bz. 


(31) 


Reflection from a semi-infinite system 

We consider an interface between two media, say vacuum and TWS, and match the incoming, the reflected and transmitted 
waves as follows 


Er = {l + yK)-\l-^)Ei„, 
E,=2{l + yP)-^Ein, 


where we have defined 





n CO £2 J 


(32) 

(33) 


(34) 


For vacuum (medium 1) Ci = Mi = 1 and for Weyl semimetal (medium 2) we assume £2 = £ and M 2 = 1 in following. 
Assuming incoming wave is polarized along x axis as E;„ = Eqx, we obtain 



The reflection and transmission matrices are given by 

1-n+n- /(«+-«_) \ 2 + n+ + n_ i{n+-n^) \ 

I i(n_ —n+) l—n+Ti- y’ I /(«_ — «+) 2+ «+ + «_ J 
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where A = 1/(1 +«+)(! +«-) and 


n± 


, 2a c 
£2 ±- b, 

K CO 


( 36 ) 


which is Eq. (7) in the main text. The reflection matrix R is used to calculate the Kerr rotation as follows. Defining (n_ — 
n+)/(l — n_n+) = The electric held on reflection becomes 

E = £o(i + . (37) 


The electric held is elliptically polarized and major axis is rotated by 9k, where 


tan20/f 


2tj sin 0 
tj2 — 1 ■ 


(38) 


Propagation of light perpendicular to node separation 

We first analyze the propagation of light in WSM assuming the direction of propagation is perpendicular to node separation 
as shown in Fig. 1(c). Although the Figure shows the incident E is polarized along y axis, to begin with we consider a more 
general case, but still keeping the x-axis as the direction of propagation. Our goal is to derive Eq. (11) in the main text. 

Using modified Maxwell equations (27) and (30), we obtain 


V X V X E = — 


£ d^E 
dt^ 


2a, (9E 

—b X 
nc dt 


(39) 


Taking the electric held to be as E(r,f) = with wavevector k = nco/c and using the identity V x V x E = VV - E — 

V^E, we obtain Eq. (11). 

Now we specialize to the case where E is polarized along y axis (case (ii) in the main text), whose geometry is shown in 
Fig. 1(c). We must impose the boundary conditions on the tangential E-flelds, + £/ = Ey, and H-flelds, — £/ — nEy = 
y^flQ/eoJy. 

We must carefully consider what contributes to the surface current density 7/. The bulk conductivity Gyy is not responsible 
for surface currents. As shown in the main text, the longitudinal current density in a WSM is not localized at the surface. In 
contrast to good conductors, a WSM has a very long penetration depth 5 ^ e'X / ne" ^ X. 

The only contribution to the surface current density is then 7/ = CJy^Ex. Here Ex{x = 0) arises from the longitudinal held in 
the WSM already described in the main text and the off-diagonal arises from the unusual surface states. To compute 
we look at the contribution of an individual surface state, labeled by the wavevector k^, and sum up all the contributions. The 
surface state labelled by k^ is localized near the x = 0 plane with a localization length*®’given by ^ (k^) = 2b/{b^ — k^). As 
one approaches the tips of the Fermi arc at kz = Xzb, the localization length diverges, and the surface state connects as it were 
with the bulk states. Thus to compute the surface conductivity we must cut off the kj’s on the scale of the wavelength of light. 
Using the fact that each individual “layer” in k-space labeled by k^ yields a 2D quantized Hall conductivity a^^{kz) = /h, 

we obtain the surface conductivity 


^ HibX). 

J^<x 2n nh 


(40) 


Uisng the surface conductivity in eq. (40) in the boundary conditions, we can calculate the reflection r = E//E'y and 
transmission t = E/jE'y coefficients. They are are found to be r = (1 — n)/(l -fn) and f = 2/(1 -fn), where the generalized 
refractive index is given by n = n -f la^bX \n{2bX)/%^n. Therefore, the t and r coefficients are functions of the node separation 
b. 
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